Compact Skyrmions, Merons and Bimerons in Thin Chiral Magnetic Films 
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A meron is a controversial topological excitation because it carries just one half of the topological charge 
unit. It is believed that it is tightly binded to another meron and cannot be observed by isolating it. We present 
a counter example, investigating the 2-dimensional nonlinear sigma model together with the Dzyaloshinskii- 
Moriya interaction, where topological excitations are merons, bimerons and skyrmions. They behave as if they 
were free particles since they are electrically neutral. A prominent feature is that the topological charge density 
is strictly confined within compact domains. We propose an analytic approach for these compact excitations, and 
construct a phase diagram. It is comprised of the helix, meron, skyrmion-crystal, skyrmion-gas and ferromagnet 
phases. It explains quite well the experimental data recently performed in chiral magnets such as MnSi and 
FeCoSi thin films, as verifies that merons are surely basic topological excitations in the system. 



Topological excitations are endlessly fascinating. They are 
constantly under investigations in all branches of physics. 
Well known examples are vortices and skyrmions QjJ] . A fan- 
tastic topological object is a meron: It cannot exist by itself 
since it carries only one half of the topological charge unit. A 
meron was originally invented as a half-instanton in the con- 
text of quark confinement in particle physics|2|]. Later it was 
introduced as a half- sky rmion in order to account for a cer- 
tain anomalous behavior in bilayer quantum Hall effects fl. 
Though the meron is a theoretically useful concept, it re- 
mains to be a long-standing problem whether it exists in re- 
ality. There is no clear evidence for its existence in spite of 
experimental endeavors |0]. 

Magnetic thin films are ideal systems to investigate and test 
various intriguing ideas on topological excitations. Indeed, 
a skyrmion crystal [4] as well as a single skyrmion [|5|] have 
been observed in chiral magnets such as MnSi and FeCoSi 
thin films. Furthermore, magnetic domains observed in ferro- 
magnets such as a TbFeCo thin film [6] are shown to be giant 
skyrmions [7] as large as ~ 1/xra. In this paper we point out 
that merons are also detectable as almost isolated objects. 

The ground state of a chiral magnet is a helical state in the 
absence of external magnetic field. The spin texture of the 
helical state has a stripe-domain structure, where the width of 
a stripe has a fixed value determined by sample parameters. A 
stripe breaks into pieces as the magnetic field increases. By 
calculating the topological charge density, we show that the 
endpoint of a broken stripe has the topological charge Qsky = 
1/2. It is natural to identify them as merons. Finite-length 
stripes are bimerons [Fig HI, among which the shortest ones 
are skyrmions. 

When we talk about skyrmions in the 2-dimensional space, 
it is implicit to assume a Belavin-Polyakov skyrmion. Its spin 
texture approaches the ground- state value only poly normally 
at large distance. On the contrary, a skyrmion must be strictly 
compact in the chiral magnet since it is embedded within a 
stripe-domain structure. We propose an analytic scheme to ex- 
plore compact skyrmions, merons and bimerons. Being elec- 
trically neutral excitations, they behave as if they were free 
particles. Based on this observation we construct a phase di- 
agram of the chiral ferromagnet. It explains quite well the 



FIG. 1 : (Color online) (a) The spin texture of a compact bimeron in a 
thin chiral magnetic film, (b) A compact bimeron is composed of two 
half-disk domains and a rectangular stripe domain. Spins are pointed 
up on the boundary and down deep inside of the bimeron. They are 
pointed forward and backward in the recutangular part, and twisting 
circularly in the half-disk part. The topological charge density is 
nonvanishing only in the half-disk parts, each of which has Q S k y = 
1/2. The half-disk part is identified as a meron. The spin texture of a 
compact skyrmion is obtained simply by removing the recutangular 
part and by patching the two half-disk domains. 



recent experimental data carried out in a FeCoSi thin filmS]. 
This fact shows that the identification of merons is justified 
not only theoretically but also experimentally. 

Our system is the 2-dimensional plane described by the 
nonlinear 0(3) sigma model Hj with easy axis anisotropy, 



and the Dzyaloshinskii-Moriya interaction (DMI), 



(1) 



Hdm = D 



J d 2 xn(x) • (V x n(x)) , (2) 



where T = (l/2)zS 2 J is the exchange energy (z denotes the 
number of nearest neighbors, S the spin per atom, J the ex- 
change constant), £ is the single-ion easy-axis-anisotropy con- 
stant, and n = (n Xj n y ,n z ) is a classical spin field of unit 
length. The DMI term breaks the chiral symmetry explicitly. 
We introduce the magnetic field h perpendicular to the plane 
with the Zeeman energy Az = SgfiB/^oh, 
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(A z /a 2 ) / d 2 xn z (x), 



(3) 
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FIG. 2: (Color online) (a) Illustration of the helical state solution 
rii(x) described by the Jacobian elliptic function $5§. Here we have 
taken k — 0.3. (b) Illustration of the spin structure with an alter- 
nating up-down and forward-backward stripe-domains. Spins are 
strictly pointed up or down on the vertical solid or dotted lines. As 
one stripe we consider the region sandwiched by two solid lines, (c) 
Illustration of a compact skyrmion described by the Jacobian elliptic 
function dTTTb . It is constructed in such a way that any cross section 
passing through the center agrees with the cross section of the stripe. 
The topological charge is strictly confined within the solid circle. 



where a is the lattice constant. 

We have emphasized previously 10] the importance of the 
magnetic dipole-dipole interaction (DDI). However, since the 
DDI constant is very small compared to the DMI constant, we 
may ignore it in determining the magnetic structure of a chiral 
magnet. Indeed, the typical size determined by the DMI is of 
the order of 40nm, where the DDI is negligible. 

We start with a study of the system in the absence of ex- 
ternal magnetic field. The ferromagnetic spin state n = 
(0, 0, ±1) is a solution of the Hamiltonian #j-dm = Hj + 
Hdm with the energy 
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However, in general, this is not the ground state. It is easy 
to prove that the Hamiltonian Hj-xm allows one-dimensional 
periodic solutions, among which the one that minimizes the 
DMI term © is given by 
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n x [x) =0, n y [x) — cn ( , k 
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n z (x) =sn , k 



(5) 



in terms of the Jacobian elliptic functions [Fig|2ta)], where k 
is an integration constant with < k 2 < 1. The periodicity of 
a(x) is 21 s with 



2i s =4<K(k; 2 ), 



(6) 



where K(n 2 ) is the complete elliptic integral of the first 
kind. The periodic state © has an alternating up-down 
and forward-backward spin-stripe structure, as illustrated in 
FigEJb). It may be called the anisotropic helical state. Note 



that n z (—£s) = 1> ti z (0) = —1, n z (is) = 1- As one stripe 
we consider the region whose width is 2i$ with the center line 
being given by n z (x) = —1. 

By substituting © into the Hamiltonian Hj.^m, the energy 
of the helical state is analytically calculable, 
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2k£K(k 2 ) 
(7) 

where L x and L y are the sample size in the x and y direc- 
tion. We determine the parameter k by minimizing £?heiix with 
respect to k. It is given by solving 2E(n 2 ) = tvirD^/T. Pro- 
vided D£ ^> T, it is solved as 
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and the energy of the helical state is 
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We compare this with that of the homogeneous state ©. We 
find that the helical state has a lower energy than that of the 
homogeneous state if > y/2T. It is interesting that the he- 
lical state does not realize in the sample when the anisotropy 
is too large. 

We switch on the external magnetic field. The Zeeman ef- 
fect enforces the increase of the up- spin region. However, it 
is impossible to increase only the width of the up-spin part 
of the stripe, which is fixed to be is- The simple way is to 
split a stripe into two stripes, since it increases up-spin region 
[FigHa)]. Let us cut one stripe at y = and then put a cap 
so that the spin field on the cross section smoothly approaches 
the up-spin value at the boundary in order to optimize the en- 
ergy [FigQIb)]. The spin texture of the cap must be given by 

n x (r, 0) — — cn I — , k 2 \ sin t 

V < 

riy (r, 6) = cn ( — , k ] cos u , 



(10) 



n z (r,0) = sn( r -^^^ 2 



for the half-disk region (r < Is, < 9 < tt) in the cylindrical 
coordinate, since it agrees with (0) at y = 0, where cos 6 = 1. 

A stripe may be broken into three stripes with one finite- 
length stripe [FigOa)]. The spin structure of a finite-length 
stripe is illustrated in FigQ] The shortest stripe is a cylindrical 
symmetric domain [FigsOa) and[2c)], whose spin texture is 
described by dTUb for the disk region (r < is, < 6 < 2tt). 

The use of a continuum approximation and of classical 
fields to represent the spins is justified as far as we analyze 
phenomena whose characteristic wavelength is much larger 
than the lattice constant. In this regime there exists the topo- 
logically conserved charge, that is the Pontryagin number, 
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FIG. 3: (Color online) (a) Illustration of compact skyrmions, merons 
and bimerons embedded in the stripe-domain structure, (b) Illustra- 
tion of the topological charge density confined within compact do- 
mains. Dotted lines show the boundaries of stripes. 



Qsky = J d 2 x p S ky(^), with the topological charge density 

PskyO) = -1^^2 £ i3 n ( X ) ' ( d i n ( X ) X d 3 n ( X )) > C 11 ) 

where z, j run over x,y with Sij being the completely anti- 
symmetric tensor. We are able to determine the topological 
charge density for various spin textures. 

First, it is trivial to see that the stripe configuration ([5]) has 
no topological density. Then, calculating it for the half-disk 
configuration (ITOl) with < 6 < tt, we find that Q& y = 1/2. 
Similarly we find Q s ^ y = 1 for the cylindrical symmetric con- 
figuration (ITOl) with < < 2tt. Clearly we can identify them 
as a meron and a skyrmion, respectively. A prominent feature 
is that the topological charge is strictly confined within a com- 
pact domain. Hence we may call them a compact skyrmion 
and so on. In general, there appear a variety of topological 
excitations [FigOa)]. We have illustrated the corresponding 
topological charge density in Fig Ob). Let us refer to this 
regime of topological excitations as the meron phase, since 
the basic excitation is a meron that appears at the endpoint of 
a stripe. 

We now estimate the energy of a topological configuration. 
The Jacobian elliptic functions are well approximated by the 
sinusoidal functions when k is not close to 1, say, k < 0.5. 
Namely, we may approximated (0) by n x (x) = 0, n y (x) = 
s'm(kx) and n z (x) = — cos(kx), where k = This is 

the well-known expression for the helical ground state in the 
isotropic system, which is the limit £ — >> oo with k being fixed. 
The spin texture of the cap (TTOt is approximated by 

n x (r, 9) = — sin(fcr) sin n y (r, 0) = sin(fcr) cos 
n z (r,0) = — cos(fcr), (12) 

for r < it /2k. We consider the isotropic system for simplicity. 

With the use of the meron configuration (fT2l) . by integrating 
the total Hamiltonian H = Hj + #dm + Hz, it is straightfor- 
ward to calculate the energy gain when a stripe is broken into 
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FIG. 4: (Color online) (a) Phase diagram. The horizontal axis is 
the temperature T in unit of D 2 //cbT, while the vertical axis is the 
Zeeman energy Az in unit of D 2 /V. (b) Average topological charge 
density in various phases. 



two stripes, 
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(13) 

We may call it the creation energy of a meron pair. The en- 
ergy gain when a skyrmion emerges in a stripe is just twice 

of AE mein 



It follows from O that AE n 



A z > A 



Helix-SkX 



with 



a Helix-SkX 



0.27D 2 /T. 
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(14) 



Since the energy gain is negative, all stripes are spontaneously 
broken into a maximum number of skyrmions for Az > 
A^ ellx_SkX , which would lead to the formation of a skyrmion 
crystal (SkX). Namely, A^ ellx_SkX is the phase-transition point 
between the helix and SkX phases. It is concluded that there 
exists no meron phase at zero temperature. 

The SkX has been discussed in literature JH @], though 
Belavin-Polyakov skyrmions are assumed on all lattice points 
with a certain cutoff. Note that a skyrmion with a cutoff 
does not have the correct topological charge. We can fol- 
low their arguments to study the SkX with the use of compact 
skyrmions without any problem. 

On the other hand, the ferromagnet (FM) phase appears in 
sufficiently strong external magnetic field, which has only the 
Zeeman energy, E^m = — (L x L y /a 2 ) Az. By comparing 
this with (TT3K it follows that the critical Zeeman energy is 



a SkX-FM 



0.84D 2 / r > 



(15) 



so that the FM phase appears for A z > Af^™. 

We proceed to construct the phase diagram in the plane of 
temperature and magnetic field. We have it already at zero 
temperature, where there exists only the helix, SkX and FM 
phases. The meron phase appears at finite temperature, since 
its entropy is much larger than that of the helix or SkX phase. 

In determining the boundary between the helix and meron 
phases, the basic excitation is a pair of merons in the he- 
lix phase [FigOb)]. It appears when a stripe is broken into 
two pieces with the excitation energy AE mQYOns being given 
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by (fT3l) . Let N be the maximum number of the topolog- 
ical charges (the maximum number of compact skyrmions) 
that the system can accommodate. When the topological 
charge of the system is n, the entropy is given by S = 
In N\/n\ (N — n)\. Then, the free energy at temperature T 
is given by F = £heiix + nAE mQYOns - TS. Letting TV oo, 
we obtain the formula for the free energy density, 

k T ( 1 \ ^ 
/ = ^heiix + qAE memns + ( q - - j , (16) 

where / = F/N, £heiix = Eheiix/N, and q = n/N is the 
average topological charge density (0 < q < 1). It is easy to 
minimize the free energy density / with respect to q. Since 
/ is quadratic in q, it yields two lines starting from the helix- 
SkX phase-transition point at T = 0. They determines the 
boundary between the helix-meron phases and the boundary 
between the meron-SkX phases, as illustrated in FigH] 

The SkX melts into a skyrmion gas (SkG) at higher temper- 
ature. We can make a similar argument to derive the boundary 
between the SkX and SkG phases, where the basic object is a 
skyrmion in the SkX phase with the excitation energy being 

A£ sky = -£ s + - tt) A z ) . (17) 

We reach at the same formula as (IT6l) with the replacement of 
£heiix by £ S kx, and A£ mer0 ns by AE sky . In this way we obtain 
the boundary between the SkX- SkG phases and the boundary 
between the SkG-FM phases, as illustrated in FiglH Finally, 
the boundary between the meron and SkG phases is given by 
A meron-skG = D 2 j 2V for j. > ^B 2 j by comparing their 
free energies. 

The phase diagram thus constructed is characterized by the 
topological charge density q as follows: 

{0 for AEi >\k B T (Helix, FM) 

\ - ^ for \AEi\ <\k^T (Meron, SkG) , 
1 for —AEi > \k B T (SkX) 

(18) 

where AEi stands for A^ merons or AE&y. This phase dia- 
gram [FigHl captures the essential feature of those obtained 
experimentally and by a Monte Carlo simulation [5], although 
they have made no distinction between the meron and SkG 
phases. We note that the topological charge density q is ob- 
servable by measuring the Hall conductance a xy of the topo- 
logical current l 1(X lllh . a xy oc q. 

We remark that all topological excitations have topologi- 
cal charges of the same sign in magnetic thin films. Thus, 
apparently, a stripe with Q sky = breaks into a number of 
bimerons each of which carries Qsky = 1- One may won- 
der how the topological conservation holds in these systems. 
We have pointed out previously J7|] that a single skyrmion can 
be created in a thin ferromagnet film by destroying the mag- 
netic order within a tiny spot with the use of photoirradiation. 
Indeed, the topological charge is well defined only when we 



can describe the spin system by continuous classical fields. 
It means that we can break the topological conservation by 
controlling the system at the lattice-constant scale or at suffi- 
ciently high temperature. Thus, in order to create topological 
excitations, it is necessary to cool the sample from high tem- 
perature. Once they are created, their stability is guaranteed 
topologically. This would be the basic reason why a rich vari- 
ety of stable spin textures have been observed|5|]. 

We have proposed a new concept of compact topological 
excitations together with their analytic expressions, based on 
which we have explored thin chiral magnetic films. Having 
identified merons as endpoints of stripes, we have pointed out 
that they have already been observed almost as isolated ob- 
jects experimentally. The reason why merons can be observed 
is that they bear no electric charge in magnetic films. Hence, 
the length of a stripe, which is nothing but a bimeron, is a 
zero-energy mode. This allows a meron to behave effectively 
as an isolated topological object though it carries only one half 
of the topological charge unit. It is worthwhile to search for 
such merons in other branches of physics. 
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